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Introduction 

This  paper  presents  the  results  of  numerical  simulations  corresponding  to  the 
MASCOTTE  10  bar  case,  done  by  LCSR  and  SNECMA.  These  computation  concern  mainly 
the  comparison  between  two  cases  : one  case  with  a secondary  atomization  model  developed 
at  LCSR  and  one  case  without  secondary  atomization.  In  this  paper  we  will  describe  shortly 
the  secondary  atomization  model  (other  details  are  available  in  ref.  [1,2]),  the  operating  point 
(condition  of  injection,  computational  domain,  ...),  and  finally  we  will  discussed  briefly  the 
main  results  issuing  of  these  computations. 


The  secondary  atomization  model 

The  secondary  atomization  model  is  based  on  an  experimental  study  which  concerns 
the  characterisation  of  the  secondary  break-up  process  in  terms  of  break-up  regimes, 
characteristic  times  and  secondary  distributions  (for  mores  details  see  ref.  [3]). 

The  break-up  regime  is  determine  with  the  following  correlation: 

We 
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Where  DR  and  VR  represent  respectively  the  density  en  viscosity  ratios  between  the  gaseous 
and  liquid  phases.  The  value  of  the  constant  C determine  the  transition  between  the  3 break- 
up regimes  considered  in  this  study:  0.25  for  bag  break-up,  0.7  for  transitional  break-up  and 
0.9  for  shear  break-up. 

When  a droplet  is  in  one  of  these  3 break-up  regimes,  i.e.  C greater  than  0.25,  an  evaluation  of 
the  duration  for  the  break-up  time  called  initiation  time  Tjnj  is  done  with  the  relation: 

T = 0.75  T We~0M  where  T = — — 
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Finally,  after  a duration  equal  to  Tini,  the  break-up  occurs,  and  the  initial  droplet  is 
decomposed  in  secondary  fragments  following  a distribution  in  size  and  number,  as  shown  in 
the  table  below: 


Size  of  secondary  fragments  (fraction  of  initial  droplet  diameter) 


Break-up  regime 


10% 


30% 


50% 


Residual 


Bag  break-up 


93,60% 


Transitional  Break-up 


94,50% 


Shear  Break-up 


86,20% 


Operating  point 

The  specification  for  the  MASCOTTE  10  bar  case  are  respectively  10  bar,  2.1 1,  50g/s 
and  23.7g/s  for  pressure,  mixture  ratio,  mass  flow  rate  of  liquid  oxygen  and  mass  flow  rate  of 
gaseous  hydrogen.  Liquid  oxygen  is  introduce  in  the  computational  domain  with  a Rosin 
Rammler  distribution  with  D32=82(im  at  the  surface  of  a liquid  core  as  shown  in  figures  1 and 
2 . 

Computational  domain  and  liquid  oxygen  injection 

The  computational  domain  consists  in  a 2D  axisymetric  geometry  (fig.  1).  The  domain 
is  400  mm  long  and  28  mm  for  radius.  The  nozzle  is  not  modelled.  Droplets  are  introduced  in 
the  computational  domain  with  a constant  normal  velocity  (fig.2). 

Overview  of  the  THESEE  code  and  condition  of  simulation 

The  code  used  for  these  computations  is  the  THESEE  code  operated  by  SNECMA. 
This  solver  can  operate  with  2D  or  3D  configurations,  multiphase,  multispecies  turbulent 
reactive  flows. 

The  computations  are  performed  with  the  following  properties: 

- Compressible  reactive  flow 

- K-e  turbulence  model 

- EBU  combustion  model 

- Sirignano-Delplanque  vaporization  model 

- Variable  thermodynamic  properties  for  the  liquid  oxygen 

- Ideal  gas  law  and  variable  thermodynamic  for  hydrogen 

In  the  computational  domain,  four  species  (3  gaseous  and  1 liquid)  are  considered  : 
gaseous  hydrogen  (GHi),  gaseous  oxygen  (GO2)  obtain  by  vaporization  of  liquid  oxygen 
(LOX)  and  gaseous  water  (H2O).  A single  chemical  reaction  is  used  to  obtain  H2O  from 
gaseous  hydrogen  and  gaseous  oxygen: 

H,  +~0;  — » Products  H,0  + Dissociation) 

The  Cv  for  H2O  is  corrected  to  take  into  account  the  dissociation  processes  of  water  at  high 
temperature.  All  the  properties  for  liquid  and  gaseous  phases  are  function  of  temperature. 

Numerical  results 

The  results  presented  in  this  part  are  essentially  a comparison  between  the  case  with 
the  secondary  atomization  model  activated  and  the  case  without  the  secondary  atomization 
model.  The  total  CPU  time  to  obtain  all  these  results  is  about  300  hours 

Mass  fraction  contours 

On  figures  3 and  4 are  represented  the  mass  fraction  for  the  3 gaseous  species  present 
in  the  computational  domain.  With  these  figures  is  it  easy  to  see  the  apparition  of  gaseous 
oxygen  near  of  the  liquid  injection  zone  due  to  the  vaporization  and  the  rapid  disappearing  of 
the  gaseous  oxygen  and  appearing  of  gaseous  water  due  to  the  combustion  process. 

The  reaction  rate 

On  figures  5 and  6 the  reaction  rate  is  represented  for  the  two  cases  of  computation.  In 
these  figures  it  is  clear  that  the  flame  form  is  directly  influenced  by  the  secondary  atomization 
process.  The  thickness  of  the  flame  is  higher  in  the  case  of  computation  without  secondary 
atomization  model  because  of  the  bigger  droplet  and  therefore  because  of  the  difference  of  the 


vaporization  times  in  these  two  cases  of  computation.  The  flame  is  positioning  nearer  of  the 
liquid  injection  zone  in  the  case  of  computation  with  the  secondary  atomization  model 
activated. 

Radial  profiles  of  mean  temperature 

On  figures  7 and  8 radial  profiles  of  mean  temperature  are  represented  for  several  axial 
locations.  For  the  two  cases  of  computation,  the  maximal  temperature  is  obtain  for  the 
X/D1=10  location,  with  a maximal  temperature  of  about  3000K. 

Axial  profiles  of  mean  temperature 

On  figure  9 and  10,  axial  profiles  of  mean  temperature  are  represented  for  4 radial 
locations.  The  maximal  temperature  is  obtain  for  the  Y/D1=2  location,  with  a maximal 
temperature  of  about  3200K.  the  oscillation  of  mean  temperature  for  computation  with  the 
secondary  atomization  model  at  the  X/dl=3  and  X/D1=4  locations  are  due  to  a divergence  of 
computation  during  the  simulation. 

The  field  of  temperature 

The  figures  1 1 and  12  are  representation  of  the  field  of  temperature  for  the  two  cases 
of  computation.  With  these  two  figures  it  is  not  easy  to  characterise  the  difference  between 
the  two  cases  of  computation. 

Droplet  location 

The  figures  13  and  14  are  representation  of  superimposition  of  mean  velocity  field  and 
droplets  in  the  two  cases  of  computation.  The  first  remark  with  these  two  figure  is  the 
difference  of  size  for  the  droplets  in  the  two  cases  of  computation  and  the  efficiency  of  the 
secondary  atomization  process.  Secondly,  the  length  of  the  “spray”  is  higher  in  the  case  of 
computation  without  the  secondary  atomization  model  because  of  the  slower  vaporization  of 
big  droplets. 

Conclusion 

A secondary  atomization  model  has  been  developed  and  used  to  simulate  the 
MASCOTTE  10  bar  case.  Results  presented  here  show  the  difference  between  two  cases  of 
computation  which  correspond  to  a calculation  with  the  break-up  model  activated  and  a case 
without  the  break-up  model. 

The  main  difference  between  these  computations  appears  on  the  flame  structure  : 
location,  length  and  thickness.  Another  difference  is  the  liquid  location  in  the  computational 
domain  : there  is  no  more  droplets  in  the  computational  domain  after  a distance  of  about  4 cm 
from  the  injection  when  the  secondary  atomization  model  is  activated  whereas  in  the  other 
case,  droplets  already  exist  up  to  10  cm. 
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Figure  1.  The  computational  domain 


Figure  2.  Injection  of  droplets 


Figure  5.  Reaction  rate  for  computation  with  secondary  atomization 


Figure  9.  Axial  profiles  of  mean  temperature  for  computation  with  secondary  atomization 


Figure  10.  Axial  profiles  of  mean  temperature  for  computation  without  secondary  atomization 
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Figure  11.  Field  of  mean  temperature  for  computation  with  secondary  atomization 


Figure  12.  Field  of  mean  temperature  for  computation  without  secondary  atomization 
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Figure  13.  Superimposition  of  mean  velocity  field  and  droplets  for  computation  with 

secondary  atomization 


Figure  14.  Superimposition  of  mean  velocity  field  and  droplets  for  computation  without 

secondary  atomization 


